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Abstract 

The characteristic feature of the ground state configuration of the 
Skyrme model description of nuclei is the absence of recognizable in- 
dividual nucleons. The ground state of the skyrmion with baryon 
number 2 is axially symmetric, and is well approximated by a sim- 
ple rational map, which represents a direct generalization of Skyrme's 
hedgehog ansatz for the nucleon. If the Lagrangian density is canon- 
ically quantized this configuration may support excitations that lie 
close and possibly below the threshold for pion decay, and therefore 
describe dibaryons. The quantum corrections stabilize these solutions, 
the mass density of which have the correct exponential fall off at large 
distances. 
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1. Introduction 



The ground state solutions to Skyrme's topological soliton model for 
baryons [|1| with baryon numbers that are larger than 1 have intriguing geo- 
metric structures with striking polyhedral symmetry p|. The simplest exam- 
ple is the system with baryon number 2, which has axial symmetry [^. The 
deuteron may be viewed as a quantized version of this ground state configu- 
ration [Q. It remains an open issue whether direct semiclassical quantization 
of these ground state configurations represent physical systems 0, ^ . 

The numerical construction of these multiskyrmion configurations is a 
demanding task |]2|, |^ . Fortunately is is possible in many cases to find simple 
rational maps, which provide remarkably accurate approximations for the 
multiskyrmion ground state configurations 0. Such rational maps may be 
viewed to represent direct formal generalizations of Skyrme's original hedge- 
hog ansatz for the the system with baryon number 1. Employment of the 
rational map approximation greatly simplifies the study of the quantized 
modes of the multiskyrmion systems [Q. In the case of B = 2 the ground 
state solution is 6-dimensional, whereas the physical solution in the deuteron 
channel corresponds to a mode in a 12-dimensional space. We here consider 
the possibility for dibaryon solutions other than the deuteron, which may be 
represented as spin-isospin excitations of the ground state solution. 

The approach adopted here is to canonically quantize the rational map 
ansatz for the ground state solution of the baryon number 2 system in repre- 
sentations of arbitrary dimension for SU{2), using the method of refs.Q, [1^, 



ri| . The corresponding states have I = J, and therefore represent dibaryon 
states other than the deuteron. By the conventional semiclassical quanti- 
zation method the lowest one of these states is very deeply bound, with 
vibrational excitations that for J = may lie below the threshold for decay 
to two nucleons and a pion 0. When canonical ab initio quantization is 
employed the ground state moves to or above the threshold for two-nucleon 
decay for all representations. In this case the representations of lowest or- 
der admit several states with / = J = 0, 1, which fall below the threshold 
for pion decay. From the phenomenological perspective the Skyrme model 
thus can accomodate narrow dibaryon states, although it does not demand 
their existence. The overwhelming experimental evidence is against the exis- 
tence of narrow dibaryon states, although some intriguing signals have been 
recently been seen in the missing mass spectrum of the reaction pd pX 
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This paper is organized into 5 sections. In section 2 the rational map 
ansatz for the baryon number 2 skyrmion is generahzed to representations of 
arbitrary dimension. This "biskyrmion" is canonically quantized in section 
3. In section 4 the resulting equations of motion are solved numerically, and 
the spectrum is obtained. Section 5 contains a concluding discussion. 

2. The axially symmetric soliton with B = 2 

The Skyrme model is a Lagrangian density for a unitary field U (r, t) that 
is described by any representation of SU{2). In an irreducible representation 
this unitary field may be expressed in terms of three unconstrained Euler 
angles a = (a^, a^, a^) as 

f/(x,t) =D^(5(x,t)). (1) 

Here the elements of the matrices are Wigner D-functions. In an arbitrary 
reducible representation unitary field can be decomposed into direct sum of 
functions. The Euler angles a then are the dynamical variables of the 
theory. 

The Skyrme model is defined by the chirally symmetric Lagrangian den- 
sity : 

C[U{^,t)] = -^TT{R,R>^} + ^TT{[R„R^]'}, (2) 
where the "right" current is defined as 

R^ = {d^U)U\ (3) 

and (the pion decay constant) and e are parameters. 

The rational map ansatz for the Skyrme field with B = 2 (biskyrmion) 
in the fundamental representation may be generalized to an arbitrary 
irreducible representation of SU{2) in the following way: 

g.(n-r)F«{.) ^ jj^^^^ ^ exp{2tn^ ■ iFnir)}, (4) 

Here -FR(r) is a scalar function ("the chiral angle") and fi is a unit vector, 
which may be defined by its circular components: 
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The boundary condition of the chiral angle Fii{r) at the orgin is -Fr(O) = 
71. In terms of Euler angles the generalized rational map ansatz may be 
expressed as 
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Given the rational map ansatz (Q), the Lagrangian density 
the following rather simple form 



(6) 

reduces to 
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The classical Lagrangian density depends on the dimension of the represen- 
tation j only through the overall scalar factor = + l)(2j + 1), which 
can be absorbed by a renormalization of the parameters . 

In contrast to hedgehog form for the B = 1 skyrmion the rational map 
Lagrangian for the skyrmion with B = 2 depends on both the polar an- 
gle "(9 and the radius r. Integration and normalization yields the following 
expression for the mass of the biskyrmion: 



2tt^ 
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(8) 
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Here the dimensionless variable r is defined as r = e/^r. Variation of tliis 
expression for the mass leads to the differential equation for the "chiral an- 
gle": 



, sin^FD(r)\ ^.^ , , sin2FB(r) 2 ^, , , 
Kir) { 1 + 4^i^ j + 2F'^{r) ^ + ^F'^{r) 



sin2FR(r) f8 \ sin^ Fji{r 



At large distances f ^ oo this equation has the asymptotic form: 

nir) + lF'^{f) - y,Fn{f) = 0. (10) 

This equation differs from the corresponding asymptotic equation for the 
chiral angle of the B=l hedhehog solution only in the value of the coefficient 
in the numerator of the last term on the l.h.s. (4 instead of 2). The solution 
of this asymptotic equation (|1^) is 

Fn{f)=Cf-'^. (11) 

The fall off rate is somewhat larger here than in the case of B=l, as the 
power of f in (12) is -2.56 whereas it in the case B=l is -2. 



3. Quantization of the biskyrmion 

We shall employ collective rotational coordinates to separate the variables 
which depend on the time and spatial coordinates: 

t/(r,q(t)) = A(q(t))f/o(r)At(q(t))). (12) 

We consider the Skyrme model quantum mechanically ah initio and thus treat 
the generalized (collective) coordinates q(t) and the corresponding velocities 
q as dynamical variables, which satisfy the commutation relations 

r,/] = -iiir^(q). (13) 

The functions _r/"^ depend on the generalized coordinates q, the explicit 
form of which is determined by the canonical commutation relations below 
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| T3| . After substitution of ( |T^ ) into the Lagrangian density (|^) the depen- 
dence of Lagrangian on the generahzed velocities may be expressed as 



L(q, q, Fji)= ^ / C{n, q(t), q(t), FR{r)y sinddrdMip 



Here we have used the notation 

= 5^(-)'"C;(-)(q)«a^(F^)c2-'")(q) 

m 

2(ao cos^ + oi sin^ g^) 2aocosg^ 

2ai 

2ao cos 2ao 



(14) 



(15) 



The functions of dynamical variables Ca'^\q) are defined in |TI|. 

Because of the axial symmetry of the rational map configuration there 
are only two different moments of inertia, which may be defined as: 



Oo — RttoiFR) — 

-4- r drr' sin' ({12 - 37r)(l + Fg) + . , 

Jit Jo V ^ 

ai = ijai(FR) = Ra_i{Fji) 

7T 



(16) 



3^263/, 



J^" drf' sin^ Fn (^37r(l + F'^) + 16^^^^^ . (17) 



The generalized momentum operators, which are canonically conjugate 
to q, are defined as 

dL 1 

P°^ = ^= 2 R3c.M)}- (18) 
These operators satisfy the canonical commutation relations 

[p„,g^] = -i5„;3, (19) 

from which it follows that the explicit form for the matrix i?,/"^(q, Fr) in eq. 
(H) is 
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(20) 



The group parameter manifold of SU{3) is the hypersphere S^. It is 
convenient to introduce the following angular momentum operators 0] on 
the sphere: 



J" = ^ {p/^,qt)(q)} = {/^c'^'Hq} , (21) 

The components of this operator satisfy the standard commutation relations. 
Note that in eq. (^) there is no summation over the index b. 

By some lengthy manipulation the Lagrangian (0) brought into the fol- 
lowing explicit form 



L (q, q, Fr) = -Mo(F«) - AM,(Ffi) + i 



ai 



ao ai y U 



■ (22) 



Here AMjIFr) represents the quantum mass correction, which may be writ- 
ten as 



AM,(F«) = ^AMo + ^AMoi + ^AM^. 



(23) 



Here the three terms on the r.h.s are quantum corrections that are due to 
the different moments of inertia, which may be written in the form 



AMn 



The dimensionless moments of inertia here are defined as 



AMn 



^AMn 



AM, = -^AM,. (24) 
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AMol = - v^TT t/ff^l 4^ gi^2 ^ ^ - 1) (2j + 3) sin^ Fr 

40 + 7r(4j(j + l)-13) 
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75 



(26) 
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(27) 



The normalized eigenstates with fixed spin and isospin i of the corre- 
sponding Hamilton operator are 



m, m 



47r 



^™.m'(q)|o). 



(28) 



The eigenvalues that correspond to these states give the masses of the quan- 
tum biskyrmion states as: 



Mr = Mo{Fr) + AM,(F^) + i + 1) + ( 



(29) 



Here rrit is the third component of the isospin. The chiral angle is deter- 
mined by the solution of the integrodifferential equation that is obtained by 



minimization of this expression (p9|) for the biskyrmion mass: 
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Here the mass parameters AM, and the moments of inertia are integrals of 
the chiral angle that is determined by solution. 
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At large distances the equation (^) reduces to the asymptotic form 



f^F'^ + 2f F^j - (4 + m^f 2) Fr = 0. 



(31) 



Here the quantity m is defined as 
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(32) 



This mass parameter describes the behavior of the chiral angle at infinity: 



(33) 



The quantity m = efT^fh represents an effective pion mass, which governs the 
asymptotic fall off exp(— 2mr) of biskyrmion mass density. In the case of 
the B = 1 skyrmion the corresponding asymptotic fall off is exp(— mr) and 
represents the Yukawa form of the pion cloud around the nucleon. 



4. Numerical results 

We have solved numerically the integrodifferential equation ( pOf ) for the 
chiral angle of the rational map ansatz for the biskyrmion (P , which provides 
a good approximation to the ground state solution for the B = 2 skyrmion. 
The resulting biskyrmion mass values are given in Table 1. In the numerical 
calculation we employed the same values for the two parameters of the the 
model, fn and e, as were obtained in ref. [0 by fitting the empirical values 



for mass and the the isoscalar radius of the nucleon. We also solved the eq. 
(|^) for the classical case using the same parameter values. The calculated 
values for the nucleon mass in the classical case are 978 MeV, 1028 MeV and 
1090 MeV respectively for the three sets of parameter values used below to 
reproduce the empirical value 939 MeV with canonical quantization. 

In the classical treatment the rational map ansatz leads to a deeply bound 
biskyrmion solution with I = J = that is stable against decay to two 
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nucleons. This is indicated the negativity of the parameter AM = M — 2M^, 
where M is the biskyrmion mass and is the calculated mass of the 
nucleon. In the canonically quantized case this state moves up to or above the 
threshold for nucleon decay. In the fundamental representation the I = J = 
state is marginally stable against decay to two nucleons with the present 
choice of parameter values. This is not the deuteron state, which has spin 1 
and isospin 0. In view of the very small energy by which this solution falls 
below the sum of two nucleon masses (- 5 MeV) and the absence of any such 
state in all other representations of SU{2) we believe interpret this result to 
be an accidental consequence of the approximate character of the rational 
map ansatz, and thus that, within the margin of error, there is no bound 
state. 

The i = I = J = state is found to lie below the threshold for decay into 
two nucleons and a pion in the canonically quantized case in both the fun- 
damental representation and the three dimensional representation [j = 1). 
There are three states with ^ = / = J = lin the j = 1 representation, 
which lie below the threshold of pion decay of the dibaryon. With the ra- 
tional map ansatz the energies of the state with / = J = is found to be 
roughly 1950 MeV and the energies of the states with I = J = 1 and isospin 
projection mt = at 2000 MeV and with rrit = ±1 at 2010 MeV respectively. 
In the 4 dimensional representation j = 3/2 all these states are found to lie 
well above the threshold for pion decay, when the same parameter values are 
used. Should the recent empirical indications [|r^ for three "supernarrow" 
dibaryons at 1904 MeV, 1926 MeV and 1942 MeV be confirmed, they could 
thus be accomodated within the Skyrme model framework, as long as their 
spin and isospin equal 0, 1, 1 respectively.. The 7 decay pattern of those 
states, do however suggest that they all have isospin 1 and = 1"^ [0, |16|. 

Typical contours of constant classical and quantum mass density of the 
dibaryon are plotted in Fig.l for the case j = 1/2 and £ = 0. The correspond- 
ing typical contours of constant baryon number density are shown in Fig.2. 
Comparison of the classical and quantum solutions in Figs.l and 2 show that 
the quantum soliton is concentrated is more compact as the densities fall off 
exponentially at infinity in for the latter. 



5. Discussion 
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Table 1: Dibaryon parameter dependence on SU{2) group representation 
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The present work represents an exploratory calculation of the possible 
quantum excitations that the the axisymmetric ground state configuration of 
the B = 2 skyrmion may support, once self consistent canonical quantization 
is imposed. We have here considered the possible states with I = J, which 
represent dibaryons other than the deuteron, which has been considered in 
ref. i. 

The calculation is based on the rational map ansatz of ref. 0, which 
provides a good approximation to the actual ground state solution. The 
calculation, within the margin of error that is associated with the rational 
map approximation does not yield states that would be stable against decay 
to two nucleons. One dibaryon state with I = J = that is stable against 
decay into two nucleons and a pion appears in both the fundamental two 
dimensional as well as in the the three dimensional representation. The 
three dimensional representation also may accomodate two such dibaryon 
states with I = J = 1. It is of course intriguing that there are some recent 
empirical indications for such states |13[, and that for those 1 = 1 would be 
favored by the photon decay pattern |T6|. 

The possibility for a dibaryon state with 1 = has drawn considerable 
experimental and theoretical interest over the past 10 years, but no confirmed 
evidence for such has yet been found [0. Should in the end no such dibaryon 
state be found, it would put constraints on the choice of representation to 
be used with the rational map ansatz in the canonically quantized case, and 
require employment of dimension greater greater than or equal to 4, as in 



12 



such the dibaryon states he above the threshold of pion decay. 
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Figure Captions 

Fig.l Typical contours of constant classical (dashed line) and quantum 
(solid line) mass densities of the biskyrmion j — ^, i — case. 

Fig. 2 Typical contours of constant classical (dashed hne) and quantum 
(solid line) baryon number densities of the biskyrmion j — £ — case. 
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